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Abstract 

The dispersion relation of the single hole in Cu02 plane is calculated in 
the self-consistent Born approximation for the thee-band Hamiltonian. We 
find that direct oxygen-oxygen hopping removes the strong anisotropy of the 
hole spectrum around the band minima. Our results compare well with re- 
cent photoemission measurements of single hole dispersion relation in the 
SriCuOiCh. 
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Typeset using REVT^ 



Recent angle resolved photoemission measurements on different copper-oxide materials 

have found a number of interesting features in the quasiparticle dispersion. The most 
interesting is the finding of an almost dispersionless (fiat) regions around (tt, 0) and (0, vr) 
points of the Brilloin zone 0. These fiat bands were observed in doped Bi2212 and YBaCuO 
compounds. Recently, APRES experiment was performed on Sr2Cu02Cl2 0, which is an 
insulating layered cooper oxide. The experiment has the virtue of measuring a single hole 
dispersion relation in a quantum antiferromagnet. It was found that the overall bandwidth 
(~ 2.2J) and the position of the minimum of the hole band at (7r/2,7r/2) agree with the 
t — J model predictions These observations leave no doubts that we are dealing with a 
hole coupled strongly to the antiferromagnetically interacting localized spins. 

At the same time some details of the experimentally found spectrum disagree with these 
of the t — J model. First of all, dispersion along the magnetic Brilloin zone boundary 
((0,7r) (vr, 0) line) is as strong as along (0,0) (tt, vr) line. Second, it turns out that 
dispersion along (0, 0) —>■ (tt, 0) line is very weak, i.e. the band is almost fiat, and the hole's 
energies at (0, 0) and (vr, 0) points were found to coincide within the experimental error [Q. 

To understand these experimental results we study the motion of a single hole in the 
three-band model of Cu02 plane. The Hamiltonian of the model is |^ 

H = J ^ S„i Sm+2g + 2(ta + tfe) ^ ^ Sm 0.m,g,a^^y^rn,g' ,a' 
"^'9 {m,g,g') a.u' 

+ ^ 2 " ^ ^ '^m,9,a«m,9',<7 - \tpp\ ^ + H -C) (l) 

{m^g^g') ^ {m,m' ,g,g') ^ 

Fermion operators a are defined on the oxygen sites, s^y = 2^0,0' , Sm is a localized copper 
S = 1/2 spin, J2{m,g,g') represents independent sum over the four oxygen sites surrounding 
cooper site m, while J2(m,m',g,g') means sum over nearest neighbors on the oxygen lattice. 
Effective hoppings ta and tb are due to the hybridization tpd between O and Cu orbitals in 
the original model f^, whereas tpp is an independent microscopic parameter. 
With the help of the canonical transformation 

bk,a = — ((1 + e'''^) afc,^,, + (1 + e'''y) a^^j,,,) 

Pk 
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dk,a = — fe'''"/^ cos ky/2 ak,x,a - e'^^'"^ cos K/2 ak^y,^ 



(2) 



am.g^.^e*'"'™, pfc = 2^1 + 7fc, 7fc = -(cos + cos ky) 



The hole's part of the Hamiltonian may be rewritten as 



N 



(3) 



k,(T 



-Mtppl (/i(^)(^fc,A,<x - dUk,a) + f2{k)idl^h,^ + blJk, 

k,a 

where we introduced short-hand notations for 

(1 + cos kx) (1 + cos ky) 



f2{k) 



pI 



(cos ky — cos kx) cos k^/l cos ky/2 



Pk 



Note that spin-fermion interaction term couples Cu spins only to the symmetric combination 
of oxygen fermion operators, i.e. to the 6-operators. Another, rf-band, appears only in direct 
— hopping. 

We will treat the copper spins within the spin wave approximation 



^± = S(l±e'«'^-)«m + 5(lTe 



a: 



(4) 



Magnon operators are defined in the whole original Brilloin zone, and go = ^r). The spin 
Hamiltonian is diagonalized then by the standard Bogoliubov transformation 



(5) 
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1 + i^a 



2<jJa 



1 



-sgn{'y, 



^ u;g = AJSJl-^'^ 



It is convenient, following to introduce staggered fermion operators as well 

^ s 

dta = IY1 (4« + 4f^s 4+go, J ' ^' ^ = ±^ 



(6) 



After somewhat lengthly algebra we find that H = Hq + Hint, where the free part is 



q k,cr,a' 



(7) 



Here 



ta+tb 



PkPk+qo is a free dispersion due to the a^S^ coupling in the Hamiltonian 



(3). Note that free dispersion term is absent in the t — J model. In the new basis direct 
oxygen-oxygen hopping looks a bit complicated 

Hpp = -4:\tpp\ Y i4,aCk,a - 4,aZk,a) 

k o-=+~ 

+h{k + go) ((4,+ - 4,-)(cfc,+ - Cfc,-) - (4+ - 4,~){zk,+ - Zk,-)) (8) 
+ [f2{k){zl+ + zl_){ck,+ + Cfc,_) + f2{k + go)(4+ - 4-)(cfc,+ - c,,_) + H.c] ), 

where + (— ) subscript denotes +\ spin projection. 

The advantage of this representation is that interaction part of the Hamiltonian has a 
simple form 



Hint = i Y U„,^'ik, q) {ct-g,aCk,a'Pq + ctyCk-q^^f^c, 

For completeness we give an explicit form of the interaction vertex U here 



ta + tb 



UqPk,^Pk-q,+ + VqPk,+Pk-q- UgPk, + Pk-q,+ + VqPk-Pk-q- 
y UqPk-Pk-q- + VqPk, + Pk-q,+ UqPk, + Pk-q- + VqPk-Pk-q,+ J 



(9) 



(10) 



Pp,- — Pp Pp+go' Pp,+ — Pp "I" Pp+qo- 



To calculate the hole's Green's function 



G{k,u) = [G^\k,uj)-J:{k,uj)]-\ 
Go{k,uj) = {uj-Hor' 

we use the self-consistent Born approximation 



(11) 



^a,cT'ik,uj) = —J2J2^'^Ak,Q)Guyik - q,uj - ujq)U^>^^>{k,q) 
u,u' q 



(12) 



In the presence of nonzero oxygen-oxygen hopping tpp, G{k, c<j) is 4 x 4 matrix, but the inter- 
action vertices U^^^{k,q) are 2x2 matrices, because there is no couphng between magnons 
and z-fermions. 

The self-consistent Born approximation (SCBA) consists in neglecting vertex corrections 
to the self-energy. It was found to work well for the t — J model [Q , and was applied to the 
three-band model in Refs. P,^. It predicts the band minima at (7r/2,7r/2) and bandwidth 
proportional to J as does the t — J model. The energy dispersion, in the absence of 
tpp, is very similar in the two models. The main difference is that hole has a spin 1/2. 
For momenta close to (7r/2,7r/2) the hole's spin is antialigned with the Copper spin, thus 
forming a singletlike excitation p|. However, away from minima a strong mixture of the 



triplet state appears p,10 



Omitting Hint for a moment we see that at tpp = the spectrum has a degenerate 
minima along the 7g = curve. Turning on tpp lifts the degeneracy, and the minima moves 



to (7r/2,7r/2) ||T0[. At the same time the energy at k = (0,0) coincides with that at (vr, 0). 
It is worth noting that the approximate relation E{tt/2,0) = 0.5E{tt/2,7i/2) holds. An 
analogous one-hole dispersion relation was obtained in |TT| by doing small-U perturbation 
theory of the one-band Hubbard model. 



Consider now the effect of interaction. Equation (|T2|) was solved numerically by the 
simplest iteration procedure. The Green's function at energy uj on the right hand side of 
(|1^) was used to calculate the self-energy at the next energy slice, u + Au. We choose Au to 
be 0.01 in units of ta + tb, and introduce the damping constant 6 = 0.01 to achieve numerical 
stability. To check for finite size effects, we used an integration with 17 points for the 'best 
fit' (see below), and found no practical difference with the 13-points integration. To check 
the validity of the procedure we used it to solve the self-consistent equation for the t — J 
model, and found good agreement with existing results |Q. 

Our 'best fit' to the experimental data is shown on Fig.l. Parameters of the fit are the 
following: ta = 0.19, ti, = 0.5, tpp = 0.76, J = 0.125 in eV. It is quite easy to achieve 
a reasonable agreement with experiment along (0,0) (7r/2,7r/2) (^r, 0) directions. 



The main problem is to fit experiment around the (7r/2,0) point. Our search for the 'best 
fit' mainly consisted in minimizing the deviation between Eq.([T^) and experimental data 
in the region of momenta between (tt, 0) and (7r/2,0). We found that spectrum E{k) is 
not affected by small variations of ta and tb while keeping their sum fixed. For example, 
ta = 0.07, tb = 0.62, tpp = 0.76 fits the experimental E{k) as well as the 'best fit'. We 
found that even better agreement along ky = line can be reached by choosing unrealistical 
values of parameters: ta = 0.31, tb = 0.94, tpp = 1.37. But the bandwidth gets higher for 
these values (Fig. 2). 

The spectral density A{k,u!) = ^Im{TrG{k,uj)) for the 'best fit' is shown on Fig.3. 
Quasiparticle peaks are observed at all values of k. At k = (0, 0) the weight of the quasipar- 
ticle pole is smallest, most of the weight is concentrated at the very strong peak a.t u = 
0. To see that lowest peaks of A{k,u) indeed represent quasiparticles we analyzed the 
(5-dependence of peaks at k = (7r/2, 7r/2), (vr, 0) and (vr/2, 0). In all cases peak fits well by a 
Lorentzian ^ ^^^^2^g2 0- As expected, the residue Z{k) is strongest at hole's minima, 
Z(7r/2,7r/2) ~ 0.25, whereas Z{7i/2,0) is 4 times smaller. It is even smaller then that at 
the (vr, 0) point ( ^^z{n'o)'^^ — ^) simply because quasiparticle's energy at this momentum is 
higher then that at the k = (7r/2,0). It is worth noticing that our value of Z(7r/2, 7r/2) 
is quite close to that in the t-J model Smallness of the ratio mdicates that 

experimental determination of the hole's energy around A; = (7r/2,0) can be difficult. 

Recently Nazarenko et al. reported theoretical analysis of the same problem. They 



showed numerically that the spectrum isotropy around the band minima can be explained 
within the single-band t — J model with the small value of intrasublattice hopping t' ~ J. 
However 'flatness' around (7r/2,0) was not explained. Also they (see also |13[) performed 
calculations of the E{k) based on the variational wavefunction for the three-band model, and 
reached an impressive agreement with experimental results at all values of k. Theoretical 
reason for such agreement is not clear to us. The only difference with our analysis is that they 
also included local Kondo — like coupling between Cu and O spins (J2 term in notations 
of Ref. |jTO|, or Kp^ in notations of 0]). Notice that such local Kondo-term has no effect 



on the spectrum at the mean-field level, because it is proportional to {Sm + Sm+2g) ■ (^m+g, 
which is zero on the Neel state. On the quantum level this term does affect structure of 
Eq.(|l2l) significantly, as it couples magnons with both c- and 2;-fermions. We have solved 
Eg . (p!2l) , reformulated to take into account J2 coupling, and found no qualitative changes in 
the E{k). 

To make a connection with existing estimates of the microscopical parameters based 
on the density-functional approach we present on Fig. 4 results of our calculations with 
parameters from Hybertsen et al. Clearly, the agreement with experiment is reasonable. 

To summarize, we found that experimental results of Wells et al. 0] can be reasonably 
well described within the three-band model with values of microscopical parameters in good 
agreement with existing estimates. The main result is that direct oxygen-oxygen hopping tpp 
is the largest parameter in the problem. This is necessary for explaining the isotropy of the 
spectrum arount the minima at (7r/2, 7r/2). At the same time our results show the presence 
of the quasiparticle pole in the Green's function at k = {tt/2, 0) also, whereas experimentally 
the possibility of no coherent excitations at this value of momentum can not be ruled out. 
Higher experimental accuracy is necessary to clarify this question. 

We thank B. Wells for providing us with experimental data. We also thank A.Chubukov, 
D.Frenkel, and K.Musaelian for discussions. We are grateful to R.Gooding and K.Vos for 
sending us their results and useful discussions. 
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FIGURES 

FIG. 1. 'Best fit' to the experimentally measured quasiparticle dispersion. 

FIG. 2. Example of fit with 'unrealistic' values of parameters. 

FIG. 3. Spectral function of hole at different momenta. Parameters correspond to the 'best 
fit', Fig.l. Lo is in units of ta + tb- 

FIG. 4. Quasiparticle dispersion with parameters based on density-functional calculations, 
Ref.[14]. 
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Fig.l of Starykh et al. 
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Fig. 2 of Starykh et al. 
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Fig. 4 of Starykh et al. 
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Fig. 3 of Starykh et al. 
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